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Role of Interlayer Electron Hopping for Spin Density Wave State 
in the Zero-Gap Organic Conductor 
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We investigate the formation of density waves in the zero-gap state (ZGS) which has been 
found in the quasi-two-dimensional organic conductor a-(BEDT-TTF)2i3 salt under the hy- 
drostatic pressure. The ZGS exhibits the cone-like dispersion for both the conduction band and 
the valence band which degenerate each other at the two-dimensional wave vectors, ±feo form- 
ing a zero gap. By using the extended Hubbard model with repulsive interaction we calculate 
the onset temperature of the spin density wave (SDW) as a function of the interlayer electron 
hopping, which by itself does not break the cone-like dispersion. It is shown that the SDW with 
wave number 2fco is induced by the combined effect of the interaction, the inter-band excitation 
across the zero-gap and the interlayer hopping. 

KEYWORDS: a-(BEDT-TTF)2l3, zero-gap state, interlayer electron hopping, spin density wave, hydro- 
static pressure 



1. Introduction 

Organic conductors, which have been investigated ex- 
tensively for many years, 1 display various states due 
to electronic correlation, for example, superconductivity, 
Mott insulator or charge ordering under varying temper- 
ature or pressure. 

Recently, it has been found that the 3/4-filled 
quasi- two-dimensional organic conductor a-(BEDT- 
TTF) 2 l3 salt under the uniaxial pressure along the 
BEDT-TTF molecule stacks (a-axis) 2 exhibits the novel 
electronic state described by the zero-gap state (ZGS). 3 
Such an exotic state appears when the Fermi surface is 
reduced to a point, and the valence and conduction bands 
degenerate at two momenta (±fen) called contact point. 
The linear dispersion around the contact point suggests 
the massless fermion. The existence of ZGS has been ver- 
ified by the first principle calculation. 4, 5 The fact that 
the temperature dependence of the resistivity becomes 
weak under the high pressure 2,6 has been analyzed by 
using the Born approximation 7 where the life time of 
the electron is inversely proportional to the energy from 
the Fermi energy due to the linear density of states. 

Although the ZGS is explained by the tight binding 
model, there are also experimental evidences 2 showing 
correlation effects in a- (BEDT-TTF) 2 I 3 salt. The stripe 
charge ordering, which exists perpendicular to the stack- 
ing axis (&-axis) exists at ambient pressure and under 
low pressure, 8 has been analyzed in terms of repulsive 
interactions by the mean field theory. 9-11 The origin of 
the superconductivity which appears with increasing the 
uniaxial pressure along a-axis, p a , has been asserted to 
come from the spin fluctuation in the presence of both 
the charge ordering and the Fermi surface. 12 At higher 
pressures, the superconductivity is suppressed and the 
system becomes the ZGS. 

The massless (or very light mass) fermion system 

*E-mail address: katashin@slab.phys.nagoya-u.ac.jp 



has been also found in graphite 13, 14 and bismuth, 15 
which display various properties, e. g. anomalous dia- 
magnetism, 13 ' 16, 17 the absence of backward scatter- 



18, 19 



and the half-integer quantum hall effect. 
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However, there are the following characteristics for 
the a-(BEDT-TTF) 2 I 3 salt 3,22 compared with those of 
graphite and bismuth. The contact point of the a- 
(BEDT-TTF) 2 I 3 salt moves in the Brillouin zone un- 
der pressures. The Fermi velocity depends on the direc- 
tion of the momentum measured from the contact point. 
The effective Hamiltonian 22 consists of the Pauli matri- 
ces a x ,a v , a z and an on the bases of the wave functions 
for the conduction band and the valence band at k (or 
— fen) while that of the graphite 13 consists of only a x 
and (j y on the basis of site representation. The basis of 
these effective Hamiltonian is the Luttinger-Kohn repre- 
sentation 23 where wave functions are represented by the 
Bloch wave functions at fen ( or — fen)- Further the inter- 
layer electron hopping 5 suggests new states in the ZGS 
of the present salt. 

The ZGS state of the gv-(BEDT-TTF) 2 1 3 salt under 
the hydrostatic pressure 24 exhibits the anomalous be- 
havior in the resistivity at low temperature. With de- 
creasing temperature, there is the rise of the resistance 
for T < 10 K. The temperature dependence of the resis- 
tivity is strongly influenced by the magnetic field, which 
is applied to the c-axis being perpendicular to the con- 
ducting plane. For H above 1 T, the magnetoresistance 
increases again at lower temperature (T < 10 K) after 
showing a constant behavior at the intermediate region 
of temperature. For graphite and bismuth, the magne- 
toresistance has only a hump. 26,27 Thus the increase of 
the resistance indicates an instability of the zero gap, 
which is the subject of the present paper. 

As a possible state expected in the ZGS, we study the 
spin density wave (SDW) state under the hydrostatic 
pressure, which can be stabilized by the interlayer elec- 
tron hopping. Applying the mean field theory to the ex- 
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Fig. 1. Structure of the conducting plane of a-(BEDT- 
TTF)2l3 where the unit cell is given by the dashed quadran- 
gle, and el, ■ ■ • , c4,pl, ■ ■ ■ ,p4 denote the respective bonds corre- 
sponding to the transfer energy of the electron hopping. 



tended Hubbard model, the transition temperature Tq 
for the SDW state is calculated where the transfer ener- 
gies for the hydrostatic pressure are estimated from data 
of the uniaxial strain and the lattice distortion under 
the hydrostatic pressure. In §2, we give formulations for 
estimating the transfer energies under hydrostatic pres- 
sure, the charge disproportionation and the linearized 
gap equation for density waves. The onset temperature 
of the SDW state is calculated in §3. Summary and dis- 
cussions are given in §4. 

2. Formulation 

2.1 model 

The crystal structure of the a-type BEDT-TTF salt 
for the conducting plane is shown in Fig. 1, where a 
unit cell consists of four BEDT-TTF molecules (A, A', 
B and C). The x,y axes correspond to a, b axes in the 
plane, respectively, while z is the axis perpendicular to 
the plane. Transfer energies, which correspond to the in- 
traplane hopping and the interplane hopping, are respec- 
tively given by 

tior.j/3 = tA(A = cl, • • • ,p4), and t z Sij6 a p , (1) 

where i,j are the indices of the unit cell, and a and (3 
correspond to A, A', B and C in a unit cell. 

The density wave is examined by employing the ex- 
tended Hubbard model, 

7i. n. ,ct 

iot 
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where a\ aa , (a =|, J.) denotes a creation operator of the 
electron and the first term is the hopping energy (eV). 
n.n. denotes the nearest neighbor sites for both intra- 
plane and interplane. The second and the third terms 
represent the on-site and nearest neighbor repulsive in- 
teractions, respectively, where (n.n.) denotes the nearest 
neighbor site in the intraplane. The intersite interaction 
Vi a -j0(= V c , V p ), which is considered only in a plane for 
the simplicity, is taken as V c for cl, • • ■ , c4 bonds and V p 
for pi, ■ ■ ■ ,pA bonds. 

2.2 transfer energy under hydrostatic pressures 

We estimate the intraplane transfer energy (Fig. 1) 
under the hydrostatic pressure by assuming the relation, 

t A = t A (0)(l + K A p a + K b A p b ), (3) 

where p a (j>b) is the pressure along a(b) direction and 
t A (0) corresponds to the energy at ambient pressure. 
Since t A is not known directly, coefficients K A and K b A 
are estimated by using the data of the uniaxial strain. 
When the uniaxial strain p a (p b ) is applied along a(b)- 
axis, the pressure p' b (p' a ) along the &(a)-axis is also added, 
where 

Pb = r lPa, 

Pa = r 2 Pb- (4) 

Such an additional pressure is needed to cancel out the 
Poisson's effect, i.e., to retain no lattice distortion along 
the &(a)-axis. Coefficients r\ and ri are estimated from 
the data of lattice parameters as r\ — 0.335 and r-i — 
0.257, which are derived in Appendix. 

Using eqs. (3) and (4), K a A and K b A are obtained as 

Ut)=i^U t)(3). <*> 

where 

T a = 1 dt A {p a ) b ^ 1 dt A (p b ) 

hA -t A {0)' d Pa ' LA ~t A {Q)' d Pb ■ {b) 

Quantities L a A and L b A denote the variation of transfer 
energies with respect to uniaxial strain p a and p b , re- 
spectively. Equation (6) is calculated from the transfer 
energies under the ambient pressure, a-axis strain at 2 
kbar and 6-axis strain at 3 kbar 3,28 (see Table I). Sub- 
stituting ri, f2, L a A and L b A into eq. (5), we obtain the 
transfer energy, t A (p) = t A (0)(l + (K% + K b A )p) for the 
hydrostatic pressure, where K A + K b A are given in Table 
I. 

2.3 charge disproportionation 

Next, we calculate the density 
(l/N)J2 i<J (al aC r a iacr)) by using the self-consistent 
Hartrec approximation, which gives the nonmagnetic 
state and the ZGS. The charge disproportionation 
corresponds to the spatial variation of (rice)- In a way 
similar to previous studies, 10-12:22 the self-consistent 
equations for the density at the respective site a are 
written as (fe = (k x , k y ,k z )) 

K)=2x:M a7 ( fc )i 2 expMfc i )/T]+l ! ( 7) 



Full Paper 



J. Phys. Soc. Jpn. 



Full Paper 



S. KATAYAMA et al. 3 



Table I. Transfer energy per unit pressure, which is normalized by <a(0). L a A and L b A are obtained in Ref. 3 based on Ref. 28. 
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cl 


c2 


c3 


c4 


pi 


p2 


p3 


p4 


a-axis strain (L A ) 


0.167 


-0.025 


0.089 


0.089 


0.011 


0.0 


0.0 


0.032 


6-axis strain (L A ) 


0.042 


0.133 


0.167 


0.167 


0.024 


0.022 


0.053 


0.032 


Hydrostatic pressure (K A + K b A ) 


0.166 


0.077 


0.194 


0.194 


0.026 


0.016 


0.039 


0.042 



(a) 



k+Q[ \k 




Fig. 2. Feynman diagrams for Aq^ and A^ in cq. (18) (a), and 
that for A^ z a in eq. (18) and A§ a in cq. (19) (b). 



where 

4 



E( £ «M fe ) + - M<W) d Pl {k) = £ 7 (fc) d ai {k), 
0=1 

(8) 

= jj E t ic:jP^~ iklri ~ ri) +2t z cosk z 5 af 3 , (9) 

(n.n.) 

U(n a ) , 1 



E a(i 



t„ = 



2 TV 



+ T7 E V i<x;j0'( n 0')- 



(10) 



(n.n.) 



AT is the number of the unit cell. For eqs. (7) and 
(8), we used the Fourier transformation expressed as 
a%aa = (1/y/N) J2k c kao£~ lk ' ri ■ We take the lattice con- 
stant as unity. The band index 7 is taken as the descend- 
ing order of £ 7 (fc), i.e., £i(fc) > &(fc) > > &{k). 
The quantity \i is the chemical potential determined by 
a condition, J2 a ( n a) = 6, due to 3/4 filling. The effect 
of t z on the charge disproportionation and \i is negligi- 
bly small, e.g., (Sn a ) ~ and Sfi <~ 10 -5 cV in the 
present calculation of £ 2 <~ 0.003. The temperature (T) is 
taken as T — > in calculating the charge disproportion- 
ation, since the T dependence of (n a ) is negligibly small 
at temperatures for the onset of SDW. The parameters 
of U, V p , V c and p are chosen to obtain the ZGS state in 
the presence of charge disproportionation. 

2.4 transition temperature for density waves 

Now, we calculate the transition temperature, Tc, cor- 
responding to the onset temperature for the density 
waves, where the operator of the spin density waves, 
(S Q Jk), S y Qa (k), S y Qa (k)) and that of the charge den- 
sity wave (CQ a (k)), are defined, respectively, as (Q = 
(qx,q y ,q z )) 

SQa(k) = ( c i + Q Q t c fcal + c fc+Qo4 C fc"T )/ 2 ' 
^Qa(^) = — i( c k+Qul Cka l ~ C fe+Qa4 Cfc "t)/2 ; 
^Qa(fc) = ( c fe+Qat Cfca T - C fe+Qai C fe«i)/2 , 

C Qa (k) = (4 +QaT c feaT + cl +Qal c kal )/2 . (11) 



In terms of eq. (11), the second and third terms of eq. 
(2) can be rewritten as 

^2 Ua ia1 a ial a i<xl a iocl 
iot 

= -ji E u(si Qa (k')s x Qa (k) + sy_ Qa (k')s y Qa (k)) 

kk'Qa 

= a7 E U(C- Qa (k')C Qa (k) - St Qa (k')S* Qa (k)) 



kk'Qa 



(12) 



and 



(n.n.),^^-' 

= F E W aft (Q)C- Qfl (k')C Qa (k), (13) 

fefe'Qa/3 

respectively. The matrix elements of V a p(Q) are given as 
V aa (Q) = 0, 
Vi 2 (Q) = K(l + e-^), 
^i 3 (Q) = ^(l + e^), 

MO) = ^.(l + e i *")» 
V 23 (Q) = ^(e^ + e i (9.+9»)) ) 

^24(Q) = ^(l + e^), 

^34(Q) = K(l + e^), (14) 

and V^/3(Q) = Vg Q (Q). Applying the mean field approx- 
imation, eq. (2) is expressed as 

Hmf ^^2^(k)c f k7a c k7r r 

+ E [ A Q a S U k ) + A Q a S Qa (k) + h.c] (15) 

kQa 

or 



Hmf = ^Z,~t{k)c{ ia c kl<J 



k-f<7 



+ E [ A Q^Q«( fe ) + A§ Q C QQ (fc) + h.c] , (16) 

fcQo 

where Aq q , (£ = x, y, z) and Aq q are the order param- 
eters of the SDW and CDW state defined as 



A Q« = -^E^-Q*( fe )> 
fe 

Aq q = ^(C/^ + 2y a/3 (Q))(C_ g/3 (fe)). (17) 
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The quantity £ 7 (fc) in eqs. (15) and (16) is the kinetic 
energy, which exhibits the cone-like dispersion and the 
zero gap located between £i(fc) and ^(fc)- The quantity 
Cfe 7<J (= ^2 a d* ai {k)ckaa) is the annihilation operator for 
the particle with the wave number k and spin a in the 
band 7. 

The transition temperature Tc for the density wave is 
calculated by using the linearized gap equations for Aq C q , 
(( = x,y, z) and Aq q , which are shown in Fig. 2. They 
are written explicitly as 



h 2 (q;3,4) 



A C A C 



= Y, UI lY (Q;a,a')A s Q i, 

0/77' 

= - E (US a p + 2V afj (Qj) 
x I 71 ,(Q;/3,a')A^ a ,, 



(18) 



(19) 



where I 77 ' (Q; a, /3) is defined as 



I 7 y(Q;a,f3) = 



1 W(£r(* + G))-/(M*)) 



TV ^ 
fc 



£ 7 (fc + Q)-£y( fc ) 



x d ai {k + Q)d* ai , (k)d* 01 (k + Q)d M ,{k) 

(20) 

and f{x){— [e x / T + 1] _1 ) is the Fermi distribution func- 
tion. The transition temperature of the density wave is 
obtained from the condition X s — 1 or X c = 1. We note 
that spin response parallel to the quantized axis (2-axis) 
is the same as that of perpendicular to the z-axis due 
to isotropic properties for both interactions and trans- 
fer energies. The summation of 7 and 7' in eqs. (18) 
and (19) takes only the conduction band (7 = 1) and 
the valence band (7 = 2) since the density wave in the 
present calculation is determined essentially by these two 
bands. Hereafter, the units of the pressure and energy 
are taken as kbar and eV, respectively. We also use the 
two-dimensional component of Q, i.e., 



Q 



q z e z 



(21) 



where q — (q x ,q y ) and e z is the unit vector perpendicular 
to the conducting plane. 

3. Spin Density Wave 

Before studying the SDW state, we describe the state 
where the interlayer transfer is absent. The ZGS under 
the hydrostatic pressure is found as follows. When only 
the kinetic energy is taken into account, i.e., without the 
Coulomb interaction, the Fermi surface exists for < 
p < 5.4 where the contact point is located below the 
Fermi energy. With increasing the pressure, the Fermi 
surface becomes small and the ZGS is obtained for p > 
5.4. The introduction of interactions suppresses the ZGS 
state. When we take a set of parameters, U = 0.4, V p = 
0.05, V c = 0.17 12 as is used in the present calculation, 
there exists the insulating state with the stripe charge 
order for < p < 2.1, the metallic state with the charge 
order for 2.1 <p< 12.5, and finally the ZGS for p> 12.5. 




Fig. 3. q-dcpcndcncc of /i 2 (<j; a, /3) ((a, (3) = (3, 4)) for U = 0.4, 
V c = 0.17, V p = 0.05 and p = 16, where the contact points 
exist on fco = (±0.970-7r, ±0.2167r). The temperature is chosen as 
T = 0.001. 
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q=2k 





V 1 

V 2 

% 22 








0.005 



0.01 



Fig. 4. T-dependence of the bare susceptibility, x^y (<?) f° r U = 
0.4, V c = 0.17, V p = 0.05 and p = 16. The cases for (7,7') = 
(1,1), (2,2) and (1,2) are plotted by solid, dashed and dotted 
line, respectively. The (2,1) component of the susceptibility is 
equal to the (1,2) component. The wave number is taken as q = 
2fc (= (-O.O6O71-, 0.432tt)). 



3. 1 density waves in the absence of t z 

We examine the property of I 7 y(q;a, (3) for t z = 0. 
The density response function of eq. (20) consists of 
the intra-band components In and I22 and the inter- 
band component 7i2(= /21) where the main contribu- 
tion comes from the inter-band one having the effect 
similar to the nesting condition. Figure 3 shows the 
q dependence of I 77 '(q) with q — (q x ,q y ) (no q z de- 
pendence due to t z = 0). Since there are maxima at 
q = ±2fe = (=F0.0607r, ±0.432vt), the incommensurate 
density wave with q = ±2k is expected from I\2{q) and 
l2i{q), i-e., the main contribution is given by an excita- 
tion from the valence band (7 = 2) to the conduction 
band (7 = 1) close to Ef around fco- 

For the clear understanding of the effect of tempera- 
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0.5 



x 



Fig. 5. z-dependence of A for the SDW (solid line) and CDW 
(dashed line) with q = 0, 2fco, where p = 16 and T = 0.001. 
The strength of the Coulomb interaction x is defined as U = 
0.4x, V c = 0.17x, V p = 0.05cc. 



ture on the 2fco density response, we calculate the sus- 
ceptibility defined as 



(<?) = - 



1 W(£r( fc + «))-/M fc )) 



£ 7 (fc + g)-£y(fc) 



(22) 



The temperature dependence of x 77 '(2fco) is shown in 
Fig. 4. With decreasing temperature, the diagonal com- 
ponent, x 77 (2feo), corresponding to the intra-band pro- 
cess, decreases and reduces to zero in the limit of zero 
temperature due to the vanishing of the density of states 
at the Fermi energy. However, the off-diagonal com- 
ponent of Xi2(2fco)(= X2i(2fco)), corresponding to the 
inter-band process rather increases with decreasing tem- 
perature. The enhancement of the off-diagonal one comes 
from the fact that the electron-hole excitation across the 
ZGS can satisfy the nesting condition. Note that, at the 
zero temperature, the off-diagonal one does not diverge 
but has a finite value. This is ascribed to the ZGS with 
the vanishing of the density of state on the Fermi energy. 
Thus it is found that the 2kg density wave is mainly de- 
termined by the inter-band pairing. 

Here we examine the effect of interaction on the den- 
sity wave by calculating A in eqs. (18) and (19). In ad- 
dition to the direct effect of interaction by U and V, 
there is also the effect through iw, which contains (n a ) 
in £ 7 . The variation of interactions are shown in Fig. 5, 
which denotes x-dependence of A with U = OAx, V c = 
0.l7x,V p = 0.05a; (0 < x < 1). The set of interactions 
for x — l 12 corresponds to the parameters, which are 
taken to explain the experimental result. 2 We note that, 
for t z — 0, such a choice of interactions (x — 1) gives 
Tc — due to A < 1 while the further increase of x{> 1) 
leads to a finite Tc with A = 1. However such a Tc may 
be reduced to zero for t z = by the two-dimensional fluc- 
tuation. Within the present mean-field treatment, the in- 
crease of Tc is noticeable in the presence of the interlayer 
transfer t, as shown later. 




(b) 



q=n 



%(k,k z +n) 




Fig. 6. Band dispersions £ 7 ' (k,k z ) (7 = 1,2) around fc = fco 
and —fco m the presence of t z for (a): q z = and for (b): q z = it, 
respectively. The Fermi surfaces are drawn by the white circle. 



In Figs. 4 and 5, A for the SDW and CDW states is 
calculated with the transfer energies at for p = 16 which 
is a reasonable pressure for the ZGS in the experiment. 24 
The fact that A oc x for small A is understood as follows. 
The main effect of interaction comes from the coefficient 
of r.h.s. of eqs. (18) and (19) while the effect on (n a ) is 
negligibly small. With increasing x in all the curve, the 
x-lincar dependence is suppressed indicating that the en- 
hancement of the charge disproportionation by interac- 
tions suppresses the density wave with q = and 2fco- 
The reason for such a behavior partly comes from the 
suppression of the density of states close to the Fermi 
surface by the increase of the charge disproportionation. 
The density response with q = 2ko is larger than that 
with q = 0. The increase for CDW as the function of x 
is due to the effect of V. 

Although the 2fco _ SDW state is the largest one among 
these four states, the magnitude of the interaction is not 
enough to obtain the finite Tc i.e. A = 1 even for x — 1 
in the present choice of parameters. 

3. 2 spin density wave in the presence of t z 

Based on the results of the previous subsection, we 
examine the SDW state, which leads to the finite Tc in 
the presence of t z . 

We note that Tc takes a maximum at q z = 7r, i.e., the 
SDW with the interlayer variation being out of phase. 
This fact has been found in the quasi-one-dimensional 
density wave system, which consists of an array of chains 
coupled with the interchain electron hopping. 30 The typi- 
cal q z dependence of the present case appears in the main 
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Fig. 7. ^-dependences of Tc for 2fco-SDW state with q z = it 
for several pressures of p = 15.6 (solid line), 16.0 (dotted line), 
16.4 (dashed line) and 16.8 (dot dashed line), where Coulomb 
interactions are chosen as U = 0.4, V c = 0.17, V p = 0.05. The 
inset figure is q z -dependences of the inter-band susceptibility 
X? 2 ( 2fc 0,<Jz) given as eq. (22) for tz = 0.0 (solid line), 0.001 
(dotted line) and 0.003 (dashed line) at p = 16, where the tem- 
perature is set as T = 0.001. 



contribution of the inter-band susceptibility x^i^Ot Qz)- 
The inset of Fig. 7 displays x° 2 (2feo7 Qz), which takes a 
maximum at q z = tt and a minimum at q z =0. The q z 
dependence of the susceptibility is essentially given by 
x" 2 (2fco,g 2 ) <x cos<7 z due to the dispersion with 2t z cosk z 
in eq. (9). Here, we note a fact that %] ) 2 (2fco J Qz) takes 
a maximum at q z = ±tt. In Fig. 6, a typical example of 
the band dispersion around the two contact points ±fco 
is shown, where £^ and & denote £ 7 (fc, k z ) (7 = 1,2) for 
k ~ fco and k ~ — ko, respectively. For q z — 0, as shown 
in Fig. 6 (a) , there is no contribution from the state with 
the wave number inside the circle of the Fermi surface, 



since electrons having ^ and £ 7 are both occupied or 
vacant, i.e. the inter-band process is absent. For q z = tt, 
on the other hand, the nesting condition is satisfied as 
shown in Fig. 6 (b), for the valence band of & and the 
electron band of £™. As for the state outside of the circle 
of the Fermi surface, it is also found that the contribu- 
tion of Fig. 6 (b) is larger than that of Fig. 6 (a). Thus, 
the inter-band process is enhanced by the appearance of 



two Fermi surfaces given by with 



and with 



k z + tt, and the SDW is optimized at q z = tt. 

In Fig. 7 the t z dependence of Tc of the 2fc -SDW state 
is shown with some choices of pressures for U — 0.4, V c = 
0.17, V p = 0.05. With increasing t z , Tc increases notice- 
ably while Tc for t z — is negligibly small. The increase 
of Tc for q z = tt by t z is ascribed to the grow of the Fermi 
surface of both valence and conduction bands, as drawn 
in Fig. 8, where the nesting condition is partly kept but 
is partly violated by t z (e.g., there are two choices of 
pairings in Fig. 6 (b)). Note that q z = tt is the relative 
momentum between £^ and £^ and then does not de- 
pend on the sign of t z . The suppression of the SDW for 
q z = comes from a fact that the region of the Fermi 




Fig. 8. Fermi surface in the three-dimensional momentum space, 
(k x ,k v ,kz), with fc z being the momentum along the interplane 
where U = 0.4, V c = 0.17, V p = 0.05, p = 16 and t x = 0.002. 
The left one (L) and the right one (L) denote the Fermi surfaces 
around —fco = (— 0.970-7T, — 0.2167r) and fco, respectively, ( for 
convenience, —fco is replaced by that of the extended zone). The 
Fermi surface is given by the valence band (v) or by the conduc- 
tion band (c), respectively, where the arrow denotes the nesting 
vector. 



surface participating in the nesting condition is always 
reduced by t z in spite of the emergence the Fermi surface. 
With increasing pressure, Tc is reduced rapidly due to 
the property of p-dependence of transfer energies, which 
increase the dip of the density of states around the Fermi 
energy under hydrostatic pressures, or the increase of the 
Fermi velocity under pressures. 

4. Summary and Discussion 

In summary, we have examined the role of interplane 
transfer energy on the SDW state in the ZGS of the 
a-(BEDT-TTF)2l3 salt under the hydrostatic pressure. 
First, we estimated the pressure dependence of transfer 
energy. The ZGS appears at pressures higher than those 
of uniaxial strain along the stacking axis. Our result is 
qualitatively consistent with the experiment. Next, we 
studied the SDW in the presence of t z and showed the 
increase of Tc with increasing t z . Since the SDW origi- 
nates from the electron-hole excitons between two Dirac 
cones with ±fco, the nesting condition is reduced due to 
the Fermi point for the two-dimensional case of t z = 0. 
Thus the appearance of the Fermi surface by t z gives rise 
to the SDW. 

Here we mention about the location of q = q max , 
which gives the maximum value of Tc- Within the nu- 
merical accuracy of the present calculation, q max coin- 
cides with 2k even for t z 7^ although the deviation 
is generally expected. Actually q max becomes slightly 
different from 2fco for the case of the Zeeman energy, 
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instead of t z (fig = 1), where 



J. Phys. Soc. Jpn. 




Qx q 0x +0.02t 



Fig. 9. qr-dependence of A Sx 'K in the presence of the Zeeman 
energy. Parameters are set as U = 0A,V c = 0.17, Vp = 
0.05, p = 16, H = 0.003 and T = 0.001, where q = 2fc = 
(-0.0607T, 0.4327T ) and q max = (-0.058tt, 0.428vr), respectively. 
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lation. 5 For the moment, there is no data of extended 
Hiickel method for interlayer hoppings, and the role of 
the three-dimensionality remains as a future problem to 
clarify the properties of ZGS in a-(BEDT-TTF) 2 I 3 salt. 
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sgn(cr) = 1(— 1) for a =\ Q). In the presence of H , X s "' v 
increases but X s * decreases. Note that the effect of H on 
T c of S x , y is nearly the same but is slightly large com- 
pared with t z . Figure 9 shows an example of the contour 
plot of A of S XiV on the plane of q = {q x ,<]y) close to 
q — 2ko(— (— 0.060-7T, 0.4327r)), where the parameters are 
chosen as U = 0.4, V c = 0.17, V p = 0.05 with p = 16 and 
Hq = 0.003. The location of the maximum is shown by 
the cross. Although the valence and conduction bands 
coincide at ±fe , q m ax(= (-0.0587T, 0.428vr)) is slightly 
different from 2fco- Since the Dirac cone is anisotropic, 
the center of the Fermi surface with an ellipse is not lo- 
cated on fco . Thus we can see the optimum wave number 
is near q = 2ko but is not exactly the same. On the other 
hand, q max for t z ^ is less effected since t z cos k z gives 
contribution corresponding to both ±Hq. 

We discuss the estimation of transfer energies and ZGS 
under the hydrostatic pressure. As shown in the previ- 
ous section, the ZGS state under the hydrostatic pressure 
emerges at higher pressure compared with ZGS under 
uniaxial pressure p a . This result is consistent with the 
experiment, since the temperature-independent behav- 
ior indicating the ZGS is observed at p a — 10 (uniaxial 
pressure) and at p = 20 (hydrostatic pressure). However 
our result for ZGS state is still under smaller pressure. 
At high pressures with p > 10, 29 it is expected that the 
distortion is suppressed due to the deviation from the 
Hooke's law. Thus the increase of tA by the pressure 
could be suppressed for that region and the ZGS is real- 
ized under higher pressure than our estimated one. 

Finally, we note on the interlayer hopping. For sim- 
plicity, t z is taken for the hopping, in which the elec- 
tron moves into a site with the same index a{— 1, ■ • • ,4) 
of the adjacent layer, and then the position of the con- 
tact point is independent of k z (see eq. (9)). However 
the hopping into other sites, which may reduce the nest- 
ing condition, is expected from the k z dependence of the 
contact point as obtained by the first principle calcu- 



Appendix: Estimation of ri and r 2 

We calculate quantities v\ and r 2 from the data of the 
lattice parameters under the uniaxial strain 28 and hydro- 
static pressure. 29 The variation of the lattice constants, 
which are defined as u a = (a(p a ,Pb) — ao)/ao for a-axis 
and lib — (b(p a ,Pb) — bo) /bo for o-axis, is written as 

f u a \ f Sl s 2 W Pa \ 

V u b J \ S3 Si J \ Pb J 

where ao = a(0, 0) and bo = 6(0,0). From eqs. (4) and 
(A-l), parameters ri, s\, s 2 , S3 and S4 are related with 
the lattice distortion per unit pressure, Ej, (j — 1, • ■ • 6). 
Actually we obtain 

s 1 +s 2 = E 1 , s 3 +s 4 = E 2 , (A-2) 

«i + s 2 ri = E 3 , s 3 + s 4 ri = E4, (A-3) 

S]T 2 + s 2 = E 5 , s 3 r 2 + s 4 = E 6 , (A-4) 

where eqs. (A-2), (A-3) and (A-4) correspond to the hy- 
drostatic pressure, a-axis strain and &-axis strain, re- 
spectively. From the experimental results of the lat- 
tice distortion, 28,29 Ei (i = l,---,6) are given as 
Ei = 0.00265 kbar" 1 , E 2 = 0.00247 kbar~\ E 3 = 
0.00326 kbar" 1 , E 4 = 0.0 kbar~\ E 5 = 0.0 kbar -1 
and E 6 = 0.00340 kbar~\ By solving eqs. (A-2), 
(A-3) and (A-4), we obtain s x = 0.00357 kbar~\ 
s 2 = -0.000917 kbar~\ s 3 = -0.00124 kbar" 1 , s 4 = 
0.00371 kbar -1 , n = 0.335 and r 2 = 0.257. 

We note that quantities n , • • • , S4 can be expressed 
in terms of Young's modulus, Y a and and Poisson's 
ratio, v a b and i>ba, which are defined as 

Y a =( P -) , Y b =( P -±) (A-5) 

\ U a) pb= o \UbJ pa= o 

Vab = —\ — \ , v ba = - ( — J • (A- 6) 

\ u bJ Pa =o \ u aj pb=0 

Since s a = 1/Y a , s 2 = -r 2 /Y a , s 3 = -ri/Y b , s 4 = 1/Y b , 
we obtain 7*1 = Vb a Yb/Y a and r 2 = v ab Y a /Yb. 
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